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Abstract. The assumption that the newly observed charged bottomonia states Zb(10610) and Zb (10650) 
are of molecular nature is confronted with the measured invariant mass distributions for the transitions of 
the T(5S) to the final states hbn + n~ and ht(2P)n + n~ . It is shown that the assumption that the Z;,(10610) 
and Zb(10650) are BB* + c.c. and B* B* bound states, respectively, with very small binding energies is 
consistent with the data. The calculation is based on a power counting for bottom meson loops, which is 
explicitly given up to two-loop in the framework of a nonrelativistic effective field theory. We also show 
that if the Z^, states are of molecular nature, then the data should not be analyzed by using a Breit-Wigner 
parametrization. 

PACS. 14.40.Rt Exotic mesons - 13.25.Gv Decays of J/ip, X, and other quarkonia 



1 Introduction 

Very recently, the Belle Collaboration reported the obser- 
vation of two charged bottomonia states .Zb(10610) and 
Zb(1065Q) in five different decay channels of the T(5S) pQ. 
Their masses and widths from averaging the measure- 
ments in various channels are Mz b = 10608.4 ± 2.0 MeV, 
r Zb = 15.6 ± 2.5 MeV, and Af z < = 10653.2 ± 1.5 MeV, 
r z , b = 14.4 ± 3.2 MeV, respectively. 

The lower Zb lies very close to the BB* threshold, 
10604 MeV, and the higher one is close to the B* B* 
threshold, 10650 MeV. The reported masses are slightly 
above the corresponding thresholds. Very soon after the 
discovery, it was proposed that the two Zb states are of 
molecular nature [2]. To be precise, the main components 
in the wave functions of the Z;,(10610) and Zf,(10650) (to 
be called Z b and Z' b ) are BB* + c.c. and B* B* , respec- 
tively. Much attention has been paid to this molecular 
interpretation. Calculations from QCD sum rules [3j, con- 
stituent quark models [3] and one-boson-exchange model 
of bottom meson potentials [5] claimed the existence of 
/ = 1, J PC = 1+- BB* and B* B* molecular states, and 
obtained masses consistent with the measured values for 
the two Z), states. In Ref. [5], the authors argue that the 
one-pion-exchange potential does not support an 5-wave 
BB* resonance state above threshold based on an effective 
field theory. However, as will be shown in this paper, the 
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experimental data in the ht,Tr + Tr~ and hb(2P)ir + Tr~ chan- 
nels are also consistent with masses slightly lower than the 
corresponding thresholds. 

As already known from many other processes and for 
a long time, coupled channels can produce peaks at their 
thresholds as a result of the unitary cut. This kind of effect 
was also noticed for the case of the Zb states in [7I8I9J . 

In fact, the existence of an isospin vector exotic state 
with J p = 1 + at the bottomonium mass region was pro- 
posed many years ago |10j . There is a long-standing puz- 
zle in the decay T(3S) — > T(1S)tttt. Before the measure- 
ments of the 7T7T invariant mass spectra of the transi- 
tions from the T(4S') to the lower-lying T states with 
the emission of two pions, two evident bumps show up 
in the 7T7T invariant mass spectrum of the decay T(3S) — > 
Y(lS)irir, see e.g. [IT]. This peculiar structure is quite dif- 
ferent from that of the transitions ip(2S) — > J/ipirir and 
T(2S) —> T(1S)tttt. Many models were proposed in or- 
der to resolve this puzzle. It was shown that by including 
the mentioned hypothetical exotic state as well as the tttt 
final state interaction, one can describe all the data of 
the T(nS) — > T{mS)Trir (n > m,n = 2,3) transitions 
well |12I13] . Later on, it was shown that the double-bump 
structure of the dipion invariant mass spectra of the tran- 
sitions T(AS) — > T(1S, 2S)irn can also be described with 
the same exotic particle [13]. However, in all these analyses 
only one state was included, while the Belle Collaboration 
reported two Zb states. Very recently, after the discovery 
of the Zb states, they were shown to play an important 
role in the helicity angular distribution of the transition 
T(5S) — > T(2S)tttt 15 . In all the above studies, effects 
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of bottom and anti-bottom meson loops which can couple 
to the exotic states were not taken into account. In view 
of the fact that these two Zb states are in the vicinities 
of the BB* and B* B* thresholds, respectively, the bot- 
tom meson loops could affect the line shapes of the Zb 
states significantly analogous to the KK loop effects on 
the /o(980) and ao(980) [16]. A systematic (re) analysis of 
all the T(nS) —> Y(mS)inr (n > m,n = 2,3,4,5) includ- 
ing both Zj, states, bottom and anti-bottom meson loops 
and 7T7T final state interaction is necessary [17] . In this 
paper, we will however focus on a simpler task, namely 
study the Zb states in the decays of T(5S) to hi,n + 7r~ 
and hb(2P)ir + iT~ which are dominated by the Zb states 
as indicated by the data [TJ , and check if the molecular as- 
sumption is consistent with the available experimental in- 
formation. Especially we demonstrate below that the data 
allow for bound state poles located below the BB* and 
B*B* thresholds, respectively. 



Fig. 1. Expansion of the two-point Green's function. Dou- 
bles lines and bubbles represent the bare propagators and self- 
energies, respectively. 



transformation properties, one can construct the Lagrangi- 
an for an 5* wave coupling of the Zb states to the bottom 
and anti-bottom mesons, 



C z = 



zl'Haa'Hb) + H.c, 



(4) 



with 2 barc is the bare coupling constant, which will be 
renormalized to the physical one z = z barc ^/Z where Z 
is the wave function renormalization constant. The three 
different charged states are collected in a 2 x 2 matrix as 



2 Lagrangians 

Because the spin-dependent interaction between a heavy 
quark and a gluon is suppressed by 1/toq, with toq the 
heavy quark mass, in the heavy quark limit, the spin of 
heavy quarks decouples. Hence, it is convenient to in- 
troduce heavy hadrons and heavy quarkonia in terms of 
spin multiplets. In the rest frame v^ = (1,0), v being 
the heavy quark velocity, one has J = T ■ a + r\b with 
T and r\b annihilating the T and hb, respectively, and 
H a = V a • a + P a , with V a and P a annihilating the vector 
and pseudoscalar heavy mesons, respectively. <j % are the 
Pauli matrices, and a is the light flavor index. Explicitly, 
one can write P a (V a ) — (i?W - , i?(*) ) for bottom mesons. 
The heavy mesons containing an anti-heavy quark are col- 
lected in H a = - V a ■ cr + P a \18\. 

The Lagrangian for the coupling of the T to the bottom 
and anti-bottom mesons can be obtained by evaluating the 
trace in the Lagrangian 



■92 



fe>- 



Cr = ^^H a a.dH a 



H.c, 



(1) 



.•h&. 



where AdB = A(dB) ~ (dA)B. The Lagrangian for the 
coupling the bottom mesons to the P-wave bottomonium 
hb reads [H] 

C hb = - gi e^ k hl^V a k +t gi hl l (V:P a ~ P a V:)+K.c. (2) 

Since the quantum numbers I G (J P ) — 1 + (1 + ) are fa- 
vored for both states by the experimental analysis [T] , the 
C parity of the neutral Zb states should be negative. This 
means under parity and charge conjugation, the fields an- 
nihilating the Zb states should transform as 



Z' 



—> z 



7iT 



(3) 



and the heavy quark spin symmetry transformation is 

given by Z % — > SZ Z S^ with S and S acting on the bot- 
tom and anti-bottom quark fields, respectively. With these 



(^Z 0i Z+* \ 

Zba = 7—j 1 7O1 

V Z -~2 Z ) ba 

The axial coupling of the pion fields to the heavy and 
anti-heavy mesons in heavy flavor chiral perturbation the- 
ory at the lowest order is given by 20 18 



r 



V2F„ 



{H\HbO ■ V<})ba) 



-jLf (H a a ■ V^ ab Hb) 



(5) 



3 Propagator of the Z b states 

The propagator of the Zb states is given by the two-point 
Green's function 



«- Gz(E ) S /^e-<«(o|T {Z ; WZ «(0)} 





' (6) 

where i,j and a,b and the indices for spin and isospin, 
respectively, and T denotes time order. Figure[T]illustratcs 
the renormalization of the Green's function to one loop 
order. The bare propagator i/[2(E — £q)] is dressed by 
the self-energy — iE. Therefore, the full propagator can 
be written as 



G Z (E) 



1 



2E-£ Q - S{E) 



(7) 



Our convention is such that the non-relativistic normaliza- 
tion differs from the relativistic one by a factor of 1/^/Mz- 
The self-energy in <i-dimcnsional space-time reads 



E(E) 



bare \ 2 



,(z barc ) 



1 



d d l 

(27r) d 1° - l 2 /(2mi) + ie 
1 



P-/°-l 2 /(2m 2 ) 



(8) 



Notice that a factor of 2 has been multiplied in this defini- 
tion in order to take into account both BB* and its charge 
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conjugated channel. The same factor appears in the B* B* 
self-energy due to a different reason: e^e^fc — 2<5" . In di- 
mensional regularization with the MS subtraction scheme, 
this integration is finite for d = 4. It corresponds to an im- 
plicit subtraction of the linear divergence which appears 
at d = 3. Taking d = 4, one has 



S MS (E) = (z*f)*JLy/Z2l&. 

an 

_ r bare \ 2 A* 



IC 



y/-2fiE9{-E)-iy/2fiE9{E) 



(9) 



where /i = mim,2/(jni + m-i) is the reduced mass, and 
6(E) = 1 for positive E and for negative E is the 
step function. The superscript MS denotes the subtrac- 
tion scheme. 

The bare energy £q is renormalized to the physical en- 
ergy, £, at the mass of the Zb- £ is connected to the mass of 
the Z\> as £ = Mz — rn\ — mi with mi 2 being the masses 
of mesons in the loop. The renormalization condition is 
such that £ is a zero of the real part of the denominator 
of the propagator. Thus, £ = £o + ReS(£), and expanding 
the real part of the self-energy S(E) around E = £ gives 



1 



Gz(E) = - 



2(E-£)[1-Re(£'(£))]-E(E) 
1 iZ 

2E-£-Z£(E) 1 



(10) 



where the wave function renormalization constant is Z = 
[1 — Re(S' (£))]" with £'(£) representing the derivative 
of S(E) with respect to E at E = £ , and 

S(E) = S(E) - Re(S(£)) - (E - £)Re(£'(£)) . 

Particularly, Re(£(£)) = Re(£'(£)) = 0. In the stan- 
dard scenario (absence of nearby thresholds; stable states) 
£(E) is dropped, however, here, due to the very close 
branch point singularity at E = 0, this function not only 
acquires an imaginary part for E > but also varies 
rapidly. It therefore needs to be kept in the propagator. 
Eq. (JTUJ) is valid for both £ > as well as £ < (but 
ill defined for £ = 0). The expression for Z follows from 
Eq. ©, 



bare \ 2 



fi 2 (z barc ) 

87T7 



-£) + lx0(£). 



(11) 



If the Z b (Z' b ) is a pure BB* (B*B*) bound state, the 
wave function renormalization constant should be since 
1 — Z measures the probability of finding a bound state 
in the physical state, see e.g. Refs. [3T]. This means the 



bare coupling 



v barc 



goes to infinity. However, the physical 



effective coupling is finite 



bare \ 2 



(z eB Y = lim Z(z baie ) 



M 



8tt 
H 2 



(12) 



with the binding momentum 7 = y/—2fj£. Eq. (|12[) co- 
incides with the one for an S wave loosely bound state 



T(55) 



(a) 



(')) 



Fig. 2. Two loop diagrams for the subprocesses T(5S) — > Zbn 
(a) and Zb — > kbit (b). Solid lines in the loops represent bottom 
and anti-bottom mesons. 



derived in Refs. [5T] taking into account the factor 2 as 
discussed below Eq. ©. 

Furthermore, the Zb states can also decay into channels 
other than the bottom and anti-bottom mesons, such as 
T(nS)ir (n — 1, 2, 3), hb(lP, 2P)n, r]i,p and so on. For the 
complete propagator we therefore need to write 



G Z (E) 



1 



iZ 



%E-£- ZS(E) + irP h v s (E)/2 



(13) 



According to the power counting in the non-relativistic 
effective field theory analyzed in details in Ref. [T5], the 
transition amplitude between the Zb states, which couple 
to the bottom mesons in an 5*- wave, and S-wave bottomo- 
nia should scale as q 2 /(M^vb), with q being the exter- 
nal momentum, and vb the velocity of the intermediate 

1/2 
bottom mesons. Because q <Si Mbv^ , this is a suppres- 
sion factor. On the other hand, the transition amplitudes 
to P-wave bottomonia hb(lP,2P) scales as q/vs- Since 
Db < 1, one would expect the decays from the Zf, to 
hb(lP, 2P)ir dominate those to Y(nS)ir. Therefore we as- 
sume that _TP h y s is saturated by the former channels. 



3.1 Power counting of two-loop diagrams 

So far we have only considered one-loop diagrams. There 
can be more loops by exchanging pions between (anti-) 
bottom mesons. Following the formalism set up in Refs. J19I 22 23 , 
we can analyze the power counting of these higher or- 
der loops. For processes with intermediate heavy meson 
loops, if the virtuality of these intermediate heavy mesons 
is not large, their three-momenta are small compared with 
their masses. Hence these heavy mesons can be dealt with 
nonrelativistically, and one can set up a power counting 
in terms of the velocity of the heavy mesons, vb- In this 
power counting, the momentum and energy of the interme- 
diate mesons scale as vb and v B , respectively, and hence 
the measure of one- loop integration scales as J d A l ~ v B . 
There are two different topologies to be distinguished. 
On the one hand there are vertex corrections — those are 
diagrams where either one or more pions are exchanged 
or where a four-i?-meson contact operator is inserted, fol- 
lowed by a two-heavy meson propagator. Since the typical 
momentum in the loop is a lot larger than the pion mass 
and the pion couples with a P-wave, the vertices for the 
exchanged pions as well as the short ranged operator pro- 



vides a factor v B , with A > 0, to the power counting of 
the diagram. In addition, the two-loop vertex correction 
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has one more integral measure as well as two more prop- 
agators, each ~ v B . Thus in total a vertex correction 
appears to be suppressed by a factor v B +x /(v B ) 2 = v B +x . 
However, one needs to pay special attention to the two- 
loop diagrams of the kind shown in Fig. [3J where, e.g., a 
pion gets produced on one heavy meson and rescatters off 
the other one before going on-shcll. These diagrams need 
to be analyzed case by case. 

Let us first analyze the power counting of the diagram 
Fig. [Ha) for Y(5S) — > Zbir. The leading order amplitude 
for the bottom meson-pion scattering formally scales as 
(E^i+E^/F 2 , see e.g. [26], with £^1,2 the energies of the 
two pions. Due to some subtle cancellation mechanisms 
also the energy of the exchanged pion gets put on-shell in 
this vertex — in analogy to what happens in the reaction 
NN — > NNir [57]. For the numerical estimates below we 
use £V = M r(5 s) - M z (t) ~ 250 MeV. There are two 

P-wave couplings in the two-loop diagrams: the coupling 
of the T(5S) to the bottom and anti-bottom mesons and 
the vertex emitting a pion inside the loop. The Zb bottom 
meson vertex is in an S'-wave. The momenta from the two 
P-wave vertices can contract with each other, and hence 
scale as v B (recall in the one-loop case, the T(5S)BB P- 
wave vertex must contract with the external momentum, 
and hence scales as q [H]). There are five propagators, and 
each of them has a contribution of order v B 2 . Therefore, 
the power counting of the two-loop diagram of Fig. dja) 
reads 

(v%) 2 v% P. „__v%E*M B (w) 



(v 2 b) 5 



16tt 2 F 2 



M f 



4 



where the factor 1/16-7T 2 appears because there is one more 
loop compared to the one-loop case, and the chiral symme- 
try breaking scale is denoted as A x = 47rP lr . We have in- 
troduced a factor of Mb to make the whole scaling dimen- 



sionless. One may estimate vb ~ \i(M — 2Mb)/Mb — 

0.15 with M = (Af T (5S) + M z )/2 and M B the average 
mass of the bottom mesons B and B* . This is to be com- 
pared to the one-loop diagram, which scale as q 2 /(M B VB) 
- see Refs. [23119] and the previous section. Numeri- 
cally, it is of similar size as or even smaller than the 
two-loop diagram given in Eq. (|T4| . Hence, the two- loop 
diagram could be more important than the one-loop dia- 
grams. However, the same mechanisms that suppress ver- 
tex corrections should also suppress three- or more-loop 
diagrams. 

The situation for the transition Z b — > hi,n is different. 
The two-loop diagram is shown in Fig. [2jb). Now there is 
only one P-wave vertex, which is the bottom meson-pion 
vertex. All the other three vertices are in an S'-wave. Due 
to the P-wave nature of the decay Z b — > hbir, the am- 
plitude must be proportional to the external momentum. 
Therefore, the only P-wave vertex should scale as q, and 
the product of the other three vertices scales, again, as 
En. Taking into account further the loop integral measure 
and the propagators, the power counting for this diagram 
reads 



(«£) 



5 \2 



Ett 



(V 2 B ) 5 16tT 2 F 2 



qM E 



MbE* 

' A 2 ' 
x 



(15) 



7T ' 
/ 
/ 



s 



T(55) ,' Z b 




hh 



Fig. 3. Decay mechanism of the process T(5S) — > Zbn — >• 
hbnn. Solid lines in the loop represent bottom and anti-bottom 



where Mb is again introduced to render the scaling di- 
mensionless. The one-loop diagrams scale like the transi- 
tions between two P-wave heavy quarkonia. According to 
Ref. [19] . the power counting is given by q/vs, which is 
numerically much larger than the scaling for the two-loop 
diagram given in Eq. (|15p . Hence, the two-loop diagrams 
can be safely neglected for the Zb — > hbir. 

In addition, counter-terms of the kind TZ^ need to be 
considered for they are needed to absorb the divergencies 
of the loop diagrams just discussed and are expected to be 
of the same importance as the loops. The corresponding 
Lagrangian to leading order of the chiral expansion reads 



Crz b , = cTZl l a dy ab + H.c. 



t*flO, 



(16) 



The analogous counter-terms for ZbhbTr can be dropped 
here for they are suppressed compared to the one loop 
diagram. Assuming the Zb and Z' b are spin partners of 
each other [2] , one can use the same coupling constant for 
them. 

A full analysis would not only require the evaluation 
of all the diagrams mentioned above but also those where 
there are no Z\' present in the processes. This will be 
studied in a later publication [17] . Here we take a more 
pragmatic point of view and simply represent the whole 
YZbir transition by the single contact term of Eq. p^]) . 
The full transition is illustrated in Fig. [3] 

One comment is in order: the Lagrangian of Eq. (1161) 

could as well mimic a compact component of the Z\ . This 
is important especially because both the two-loop and one- 
loop scaling given by Eq. (fTH) and q 2 /(M b vb) are much 
smaller than 1, and hence a small compact component 
could be more important than the bottom meson loops in 
the TZbn vertex. Thus, since we expect this contact term 
to appear at leading order, it seems as if we would not be 
able to disentangle a compact, say, tetraquark component 
from a molecular one. However, since in the molecular sce- 
nario the transition Zb — > iihb is dominated by the loop, 
the structure can indeed be tested, since for molecules 
the dynamics appears to be quite restricted. Thus, in the 
present approach the T(5S)ir vertex provides a source 

term for the Z^ , while their decays are the subject of 
this study. 



M. Cleven et al.: Bound state nature of the Z^ states 



4 Results 



The decay widths of the Zb 
taincd to be 



and Z' b into h^n are ob- 



The relevant vertices follow from these Lagrangians. The 
decay amplitudes for the two-body transitions can be found 
in Appendix [21 On one hand, the T(5S) is only about 
120 and 70 MeV above the Z b n and Z' b ir thresholds, re- 
spectively. On the other hand, it has a large width 110 ± 
13 MeV, and the experimental data were taken in an en- 
ergy range around the T(5S) mass. Therefore, when cal- 
culating its decay widths, one has to take into account 
the mass distribution of the T(5S). Its three-body decay 
width can then be calculated using 



i r-(M r +2r Y ) 2 (0,^4 I- 

r(T(5S)h- hody = - / ds { -^ / d$M\ 2 



x— Im 

7T 



S-M2. + iM r r r 



, (17) 



where M r = 10.865 GeV and T T = 0.11 GeV are the mass 
and width of the T(5S), respectively, and f d<p3 denotes 
the three-body phase space, see e.g. [53]. The function A 
contains all the physics and can be easily obtained using 
the loop amplitudes given explicitly in the appendix. Both 
positively and negatively charged Zj, 
be considered. The factor 1/W with 



and Z' b states should 



-1 



r.(A/ r +2r r ) 2 2 

H' = / ds-Im . 

(Mr-2/Y) 2 n \S - M$ + lM Y r Y 



is considered in order to normalize the spectral function 
oitheT(5S). 

Weconsider the casethat the Zt and Z' b couple only to 
the BB* + c.c. and B*B* channels, respectively. Coupled 
channel effects should be suppressed because \£\ <C Mb* — 
Mb for both Zt, states. 

The parameters of the model are the normalization 
factors N, chosen individually for the two final states, the 
physical couplings, which are products of \/Z and the bare 
couplings, for the Zb states to the relevant open bottom 
channels, z\ and z 2 (or equivalently the binding energies 

of the Z\ ), and those couplings for the hb and hb(2P), 
denoted by g\ and g[, respectively. The parameter c of 
Eq. (|16p is absorbed into the overall normalization factor. 
Both gi and g\ only appear in a product with the Zj. In 
order to reduce the number of free parameters, we assume 
g[ = g\ in the following. Note that neither of them can be 
measured directly, since the masses of hb and hb(2P) are 
below the BB* threshold. In the actual fit we will adjust 
Z\, r z = z%/z\, g\Z\ and the two normalization constants. 

Using the amplitudes of Eqs. (IA.2IA.3[) . we fit the pa- 
rameters to the invariant mass spectra of both h b n + and 
hb(2P)n + from 10.56 GeV to 10.70 GeV in the missing 
mass spectrum MM(ir). In the chosen region, there are 
14 data points for the T(5S) — > hbir + n~ and 13 for the 
T(5S) ^ h b (2P)x+%- . 



r(Z b -> h b n) = 4.8 ( ^iGeV 2 1 NF.-V 



r(z b -> hb-K) 



140( 5 iziGeV) 2 MeV, 
5.8 (^GeV 2 1 M.V 



V F^ 
169( 5 iz 2 GeV) 2 MeV. 



(18) 



Due to smaller phase space, the widths for the decays 



A') 



hb(2P)ir get smaller numerical factors. We find 



r{z b - 

nz' b - 



■h, 
> h b 



,(2P)7r) = 30( 5 ^!GeV) ; .,„ , . 
b (2P)7r) = 46<y i z 2 GeV) 2 MeV. (19) 



) 2 MeV, 



In getting the above numbers, we use F^ = 92.4 MeV. Be- 
cause the B* is below the Bit threshold, the axial coupling 
constant g cannot be directly measured. Fortunately, there 
have been quite a few theoretical determinations using dif- 
ferent methods. For a collection of these results, see [2"5] . 
Almost all the determinations fall in the range between 
0.3 and 0.7, and g = 0.5 is used here. 

If the Zb and Z' b states are S-w&ve bound states of 
the BB* and B*B* , respectively, their coupling strengths 
to the bottom and anti-bottom mesons are related to the 
binding energies. The relation has been derived in Eq. (|12p . 
The coupling constants Zi appear in Eqs. (|18[) . which enter 
the Zb propagators, as well as in the transition amplitudes 
-c.f. Eos. (lA.2IA.3jl . 

In order to fit to the data which are events collected per 
10 MeV, we integrate the invariant mass spectra for each 
bin corresponding to the measurements. This is important 
for narrow structures. The best fit results in x 2 /d.o.f. = 
54.1/22 = 2.45. 



z\ 



0.75± ;?i GeV" 1 / 2 , 



-0.39 



+0.06 
-0.07' 



g-LZi =0.40 ±0.06 GeV 



(20) 



The results from the best fit are plotted in Fig. [4] together 
with the experimental data. Using Eq. (fl"2")) the couplings 
can be converted to binding energies. Especially we find 

-O.II+H4 MeV . (21) 



£ Zb = -4.7+ 2 ; 2 MeV, £ z . 



Although the Zb and Z' b are supposedly spin partners, 
a value of r z = —0.4 is not completely unreasonable: the 
fine tuning necessary to put a bound state as close as 
0.1 MeV to a threshold is extremely sensitive to even a 
small variation in the scattering potential, driven by spin 
symmetry violations. In effect this can give significant dif- 
ferences in the binding energies and, via Eq. (TT2"j) . also in 
the coupling constants. However, a microscopic calcula- 
tion, which goes beyond the scope of this paper, would be 
necessary to check this hypothesis. 

One can obtain a better fit to the data if one would 
either release the bound state condition given as Eq. (fT2")l . 
and allowing the masses of the Zb states to float freely, or 
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Fig. 4. Comparison of the calculated invariant mass spectra of 
ht7v + and hb(2P)ir + with the measured missing mass spectra 
MM(tt). 



allow for non-resonant terms. But this is not the purpose 
of our paper — we here only want to demonstrate that the 
data are consistent with the bound state picture, which 

implies the masses of the Z\ states to be located below 
the corresponding thresholds. 

It is interesting to look at the Z b line shape or the ab- 
solute value otGz(E) for different locations of the Zb pole. 
The function \Gz(E)\ using the parameters from the best 
fit are plotted as the solid line in Fig. [5] In this case, the Zb 
is a BB* bound state with a binding energy of —4.7 MeV. 
Keeping Zl = 0.75 GeV~ 1/2 , and g lZl = 0.4 GeV" 1 fixed 
we also plot as the dashed line the line shape for the vir- 
tual state with the same value of £ Q The dotted and dot- 
dashed lines are for a resonance with a mass above the 
BB* threshold by 8 and 20 MeV, respectively. From the 
figure, one sees that the bound state produces a bump be- 



1 The curve for virtual state is obtained using Eq. (|10p but 
with E{E) = E(E) - Re(Zn(£)) - (E - £)Re(E'(£)) , where, 
in the MS scheme, Re(Zii(£)) = 
in the second Riemann sheet. 




-Ke(E(£)) is the self-energy 



E [MeV] 

Fig. 5. The absolute value of Gz(E). The solid (red) and 
dashed (blue) are for a bound state and virtual state, respec- 
tively, with the same mass, £ = M Zb -M B -M B * = -4.7 MeV. 
The dotted and dot-dashed (black) lines are for resonances with 
£ — 8 and 20 MeV, respectively. The maxima of the resonance 
and virtual state curves have been normalized to the bound 
state one. 



low the threshold, and a small cusp at the threshold, while 
the virtual state produces a prominent cusp at the thresh- 
old and no structure below. Above the threshold, the en- 
ergy dependence of the virtual and bound state curves are 
exactly the same. The bump in the bound state case re- 
flects the pole position. Hence, if we reduce the value of 
binding energy to, say about 0.1 MeV, which is the case for 
the Z' b and the A (3872), then the bump below threshold 
would be invisible, and the cusp dominates the structure. 
In this case, it is hard to distinguish between the bound 
state and virtual state scenarios. For more discussions on 
the shape of a virtual state, see e.g. |28l29j . 

One important feature of the line shapes of dynami- 
cally generated states is shown in Fig. [5J for poles slightly 
above the threshold, since the coupling to the opening 
channel is strong, the position of the peak is locked to the 
threshold, as can be seen from the dotted line. Increas- 
ing the resonance mass, the effect of the cusp is smeared 
out, and shape is approaching a normal Breit-Wigner res- 
onance — in the dot-dashed line one starts to see a bump 
above threshold developing for a mass as large as 20 MeV 
above the BB* threshold. However, even then the peak is 
still located at the threshold. We are therefore to conclude 
that a Breit-Wigner parametrization, as was used in the 
experimental analysis, should not be used when analyzing 
structures that emerge from dynamically generated states. 



5 Conclusion 

We showed that the data [T] is consistent with the assump- 
tion that the main components of the lower and higher Zb 
states are S'-wave BB* + c.c. and B* B* bound states, re- 
spectively. A small compact tetraquark component, how- 
ever, can not be excluded. 



M. Cleven et al.: Bound state nature of the Z^ states 



It is difficult to distinguish between resonance and 
bound state scenarios with the current data, however, data 
with higher resolution should allow one to distinguish the 
two cases — see, e.g., Refs. [30131] for the corresponding 
discussion for the A (3872). In a next step we will improve 
the model by inclusion of non-resonant terms as well as 
the calculation of of other decay channels [17] . 

We have demonstrated that if the Zf, states are indeed 
generated from non-perturbative BB* + c.c. and B* B* 
dynamics, the data should not be analyzed using a Breit- 
Wigner parametrization. This statement can also be re- 
versed: if a near threshold state can be described by a 
Breit-Wigner form, it is not dynamically generated, as 
this is possible only if the coupling of the resonance to 
the continuum channel is very weak. At present the data 
appears to be consistent with line shapes that result from 
dynamical states as well as genuine ones. Therefore a de- 
cision about the nature of the Z\, states will be possible 
only once data with higher resolution and statistics will 
be available. 
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A Loop function and expressions of the 
amplitudes 

The basic three-point loop function worked out using di- 
mensional regularization in d = 4 is 

/(mi,TO2,m 3 ,q) 

-i r d d i i i 





1 




1° + bi2 - 


" &23 - 


(i-q) 2 i ■ 
i — — — h it 

7713 



/^12/^23 1 



16tt 



tan 



tan 



2y/a(c — ie) 
2a + c - c' 



2^J a(c! — a — ie) 



(A.l) 



where rrii(i = 1,2,3) are the masses of the particles in 
the loop, jj,ij — mirrij / (rrii + rrij) are the reduced masses, 
b\2 = mi + m2 — M, 623 = 1712 + m>3 + q° — M with M the 
mass of the initial particle, and 



M23 



m 3 



M23„2 



a = — q , c = 2^12612, d = 2/i23^23 H q 
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